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We consider the noncommutative Abelian-Higgs theory and investigate general 
static vortex configurations including recently found exact multi-vortex solutions. 
In particular, we prove that the self-dual BPS solutions cease to exist once the 
noncommutativity scale exceeds a critical value. We then study the fluctuation 
spectra about the static configuration and show that the exact non BPS solutions 
are unstable below the critical value. We have identified the tachyonic degrees 
as well as massless moduli degrees. We then discuss the physical meaning of 
the moduli degrees and construct exact time-dependent vortex configurations 
where each vortex moves independently. We finally give the moduli description 
of the vortices and show that the matrix nature of moduli coordinates naturally 
emerges. 
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1 Introduction 



The noncommutative solitons found in the noncommutative scalar theory [0] do not even exist 
in the commutative version of the theory. This indicates that the characteristic properties 
of solitons in some noncommutative field theories may greatly differ from those of ordinary 
solitons. Of course, there are examples where the nature of noncommutative solitons and 
the corresponding ordinary solitons are quite similar to each other in the sense that the 
properties of noncommutative solitons are given by just smooth deformation governed by 
the noncommutativity scale 9. 

One such example is the U(2) BPS monopole discussed in Ref. [|2], [| |, ||. The energy 
and the charge of the BPS monopole do not depend on the noncommutativity scale. The 
effect of the noncommutativity appears as tilting of D-strings in the transverse space giving 
dipole nature of the magnetic charge distribution. It can be argued that the interactions of 
the U(2) BPS monopoles are independent of the noncommutativity scale 9 within the moduli 
space description of their dynamics]^, [|. Contrary to monopole case, the noncommutative 
scalar solitons found in |TJ are genuinely noncommutative object since they cannot exist 
in the ordinary scalar theory. As discussed in Ref. 0, the shape deformation of the scalar 
soliton is quite peculiar when moving with a constant velocity. Specifically, their deformation 
is not simply dictated by the Lorentz contraction but described by an area preserving ellipse 
exhibiting the UV/IR mixing phenomena of noncommutative field theories. 

We here pursue a similar issue on the recently found exact multi- vortex solutions || in 
the noncommutative Abelian-Higgs theory [0, [LTfl. (For soliton solutions of some other mod- 



els, see Ref. JTT], [12], [13], [14], [15], [L6|, [Tj]].) Certain apparent properties of the noncommutative 
vortices are striking even in their static properties. The multi-vortex solutions are in general 
not BPS saturated states but their energy, nevertheless, scales linearly in the number of 
vortices. This seems to imply that there are no interactions between vortices even in this 
non BPS case. We shall show that the self-dual BPS solutions exist only when 9v 2 < 1 where 
v is the vacuum expectation value of the Higgs scalar. This property is also contrasted to 
the commutative Abelian-Higgs theory where the self-dual BPS vortices exist for all vacuum 
expectation values of the scalar. There is another aspect concerning the noncommutative 
vortex solitons; the theory allows exact time-dependent solutions of vortices each of them 
moving in an arbitrary velocity from an arbitrary initial location. In view of generic com- 
plexity involved with soliton dynamics of field theory, the existence of such time dependent 
solutions is quite peculiar. 

In these respects, the systematic approach toward the understanding of the noncom- 
mutative vortex solutions seems imminent on the following issues. First the possible static 
solitonic configurations need to be mapped out including the self-dual or anti-self-dual BPS 
branches. Second the stability of the non BPS multi- vortices is a priori unclear. This issue 
can be studied by turning on general perturbations around the static solutions. In case there 
are tachyonic degrees possessing a negative mass squared, the static configurations are nec- 
essarily unstable. Any small perturbations in this direction will make the vortices collapse 
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to a stable configuration. On the other hand, when fluctuation spectra do not possess any 
tachyonic degrees, any individual vortex works as an stable solitonic object. The massless 
fluctuation is responsible for the moduli motions. Finally, one is interested in the interac- 
tions between vortices especially when they are stable. The interaction can be studied by 
adopting the scheme of the moduli space approximation. In fact, one may go beyond the 
moduli space dynamics by identifying quartic potential depending on the moduli coordinates 
in our present problem. Denoting number of vortices by m, the U(m) matrix nature of the 
moduli coordinates emerges and the dynamics turn out to be described by the matrix model 
of m DO-branes. 



In this note, we review first the exact solutions of non-BPS multi-vortices. We also 
describe the exact solutions where vortices are positioned in arbitrary locations. In Section 
3, we study other static solutions focused on the self-dual BPS branch. Our study will be 



summarized in Figure 1 where the anti-self-dual branch discussed in Ref. JTU| is also included. 
In Section 4, we study the general fluctuation spectra around the static solutions identifying 
all the tachyonic modes and massless modes. Masses of the degrees connecting the vortex 
to the vacuum can be identified by diagonalizing the kinetic and quadratic potential terms 
simultaneously. The remaining degrees will be shown to be equivalent to the fluctuation 
spectra about the vacuum of the original Abelian-Higgs system. In Section 5, we identify 
the moduli parameters appearing in the exact solutions by analyzing the translation and the 
moments (constructed with help of covariant position operator). We then construct exact 
time dependent solutions describing vortices moving in arbitrary velocities. The moduli space 
description is then worked out and the relevant metric will be shown to be flat. We then 
describe how the matrix nature of the moduli coordinates emerges. Last section comprises 
the summary of our results and concluding remarks. 



2 Exact multi-vortex solutions 



We begin by recapitulating the properties of the exact multi- vortex solutions of the noncom- 
mutative Abelian-Higgs theory found in Ref. H. The noncommutative Abelian-Higgs model 
in 2+1 dimensions is described by the Lagrangian 

L = -j 2 j d H\ F ^ * pilv + D ^ * + ^ * <^ w ) 2 ) C 1 ) 

where 



fJ,U 



D, 



iA^c 



(2) 



The *-product is defined by 



f(x)*g(x) = (e- i ^ id ^f(x)g(x')) 



(3) 
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where we take 9 to be positive without loss of generality. The theory can be equivalently 
presented by operators on the Hilbert space defined by 

[x,y] = -i6 (4) 

where the *-product between functions becomes the ordinary product between the operators. 
For given function 

f{x > y) = j£^f( k y (kxx+kyy) > (5) 

the corresponding operator can be found by the Weyl-ordered form of 

/(*,») = / ^ (6) 

One may then easily show that / d 2 x f is replaced by 2TT0trf and dif corresponds to 
— ij€ij[xj,f]. With the operator-valued fields, the action can be written as 

L = -^tr [ -F, V F^ + D^{D^ + ^ - v 2 ) 2 ] (7) 

where hats are dropped for simplicity and the derivative notation is understood as dif = 

At this point, we introduce the creation and annihilation operators by d" = -^{x + iy) 



and by c = ^h^{x — iy), which satisfy [c, c^) = 1. To represent arbitrary operators in the 
Hilbert space we shall use the occupation number basis by G = J29ki\k){l\ with the number 
operator c^c. We will further denote A = A x — iA y , d-G = (d x — id y )G = J^[c,G] and 

d + G=(d x + id v )G = -y/%[J,C\. 

The system is invariant under the gauge transformation, 

A^rfApU + irfdnU, <j/ = rf<f>, (8) 
where the gauge group element U satisfies 

UU ] = U ] U = I. (9) 
We introduce a covariant quantity K defined by 

A = -i^(c-K), (10) 

which transforms as K' = WKU under the gauge transformation in (|8]). Later it will be 
interpreted as a covariant version of position operator up to numerical coefficient. 
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The Hamiltonian can be constructed as 

H = ^tr [\{E 2 + B 2 ) + A0(A0) t + A0(A0) t + - v 2 ) 2 ] . (11) 

g Z Z 

using the time translational invariance of the system. On the gauge choice Aq = 0, the 
equations of motion read 

<t> - DiDicf) + \(0 - v 2 )<J) = , 

Ai + eijDjB = Ji= i[(p(Di4>y - A00 f ] , (12) 
with the Gauss law constraint 

DiAi = Jo = {[<$- j><fi]. (13) 

The exact multi-vortex solutions found in Ref. are given by 

K = S m cS ] m , <p = vS m , (14) 

where S m denotes the shift operator S m = zZ^Lo \n + m)(n\ (m > 0). The shift operator 
satisfies relations 

SmSm = / , SmSjn = P m = I P m , (15) 

with the projection operator P m defined by 

m— 1 

p m =J2\ a )( a \- ( 16 ) 

a=0 

The magnetic field of the solitons reads 

B= l -P m . (17) 

The flux defined by <3> = 9ti B is m on the solution. Thus the solution describes m vortices 
of the Abelian-Higgs theory characterized by topological quantity $. The energy of the 
vortices is evaluated as 

M(M) = ^(i + A^)>^VV. (18) 

When A = 1, the theory allows so called Bogomol'nyi bound as discussed in Ref. [jTIJ. In 
fact it is straightforward to verify that the energy functional can be expressed as a complete 
squared form plus a topological term by 

-wf) ?irn 2 

H = — tr \{B ± (#t _ v 2)y + 2 (D ± 0)(D ± 0)t ± t i] D l J j ± 2v 2 B] > — , (19) 
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where we omitted the kinetic terms involving E{ and D t 4>. The saturation of the bound 
occurs once the self-dual Bogomol'nyi equations, 

D +( j) = 0, B = v 2 - 00 f , (20) 

or the anti-self-dual equations 

£>_</> = 0, -B = v 2 -<f)(fi, (21) 

are satisfied. When A = 1, the bound in fll8|) agrees with the Bogomol'nyi bound that is an 
absolute energy bound for m vortex solution. Hence when v 2 = | and A = 1, the solution 
should be a BPS solution. Indeed for the specific value of 9v 2 , one can check that the 
solution satisfies the self-dual BPS equations. This BPS solution is clearly stable because 
they saturate the energy bound set by the topological quantity. 

Other obvious generalization of the static multi- vortex solution is given by]TB[ 

i m— 1 

K = S m cS j m + —= £ A B |o)(o| , = vS m , (22) 
V 20 a=0 

where A a 's are constant complex numbers. This solution has the same flux and energy as the 
solution in fll4[) . Hence we see that A a is the moduli parameters of the multi- vortices. Later 
we shall clarify the stability of the vortex solutions, which is a priori not clear because they 
are not always BPS saturated solutions. But before discussing this matter, we will study 
the BPS solutions for A = 1 and Ov 2 ^ 1 or other possible static solutions. 



3 BPS solutions of multi-vortices 



In the last section, we have derived the BPS equations of the Abelian-Higgs theory with 
A = 1. The static multi- vortex solutions are in general not BPS saturated. However they 
become self-dual BPS solutions for a special value of 9v 2 = 1. In this section we focus on 



the BPS solutions. Some analysis on the anti-BPS solutions is carried out in [10 and the 



comparison will follow at the end of this section. In terms of K the BPS equations become 

^ (l — [K, K^]j = v 2 — <j)$ , 0^-^0 = 0. (23) 

By virtue of the explicit form of c 1 the latter can be solved 

1 00 1 

0=^E -^Kn<Po)(n\, (24) 
V0n=o vn! 

where | 0o ) = V0(j)\0) is an arbitrary constant vector. Substituting this expression, the 
BPS equations are reduced to a single equation 

OO 1 

6v 2 - 1 + [K,tf] = £ -rKTlfaHMK"- (25) 

n=0 n - 
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To solve this equation we take an ansatz for K as 



K=J2fn\n){n+p\, (26) 

n=0 



where p is any positive integer. Substituting this expression into Eq. fl25|) , one can show 
the following: BPS solutions exist only when 1 > 9v 2 . Furthermore, finite energy or flux 
solutions exist only for p = 1. Specifically we obtain 



K = ^^/o(l-^ 2 )|a-l)(a | - J2 k n I n+m-l)(n + m\, 

a=l n=l 

C ^m)(0\ + ^ lll2 '" kn \n+m)(n\) , (27) 



V9 



n=l 



where ( G C, m corresponds to the flux number which is non-negative integer, and the 
sequence, k n , n — 1, 2, • • -, satisfies the recurrence relation 

Qn+i + Qn-i ~ 2g„ = — (q n - q n -i + 9v 2 ) , (28) 
n 

with q n = \k n \ 2 — n. The initial data for the recurrence relation are 

q = m(l-9v 2 ) qt = m(l - 9v 2 ) + |C| 2 - 9v 2 , (29) 
where ( is an adjustable parameter. The magnetic field is given by 

m— 1 ^ oo 

B = J2 v 2 \a)(a\ + - (Qn - ?nfi)l« + m)(n+m\ , (30) 

a=0 " n=0 

so that the flux is 

$ = m — lim q n . (31) 

Eq. (|2T| ) satisfies the BPS equations for any value of C- Choosing C = or \(\ — V9v 2 gives 
plus infinity or minus infinity flux solution respectively. Furthermore, if q n converges, the 
converging value must be zero. Thus by continuity, there exists (, < \(\ < \J 9v 2 , which 
makes q n converge to zero, and hence BPS solutions have finite and quantized energy. The 
appendix contains our proof[[ For 9v 2 = 1, the choice \(\ = 1 leads to the exact solution 
with q n = or k n = -Jn. 

We have shown that, within the ansatz taken, there are no BPS solutions possessing 
a positive flux for 9v 2 > 1. Without limiting the discussions to the specific form, one 
may prove that there are indeed no self-dual BPS solutions for 9v 2 > 1 as an analytic 
perturbation of small parameter uo around the 9v 2 = 1 BPS solution. For this purpose, we 
shall take a generic perturbation around the solution and expand it as a power series of 



§ We here like to mention that we have also numerically verified that the value of which makes the 
series to converge to zero up to a few hundred terms, approaches to a unique value for a given value of 
6v 2 G {0.1, 0.2, • • • 0.9} and for m = 1. 
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not known 
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9v 



9v ; 



a) vortices with positive flux 



b) BPS vortices with negative flux 



Figure 1: The energy of (a) the self-dual vortices and (b) the anti self-dual vortices. We also depict 
the energy of the non-BPS vortices that have a positive magnetic flux for A = 1. 



the small parameter uj = J\e\ with e = Ov 2 — 1. (The choice u = e will quickly lead to a 
contradiction.) Namely, we consider the fluctuation around the exact solution as 

<f> = v(l + (f)S m 

K = S m cSl + h, (32) 

with the expansions, 

oo oo 

y = S^V(0> h = J2v l h(i) , (33) 
i=i i=i 

and may show that there are no solutions for Ov 2 > 1. The proof is relegated to the appendix. 

One could also try an expansion with respect to a parameter uj n defined by \e\ l l' n 
for arbitrary nonnegative integers. Though a little complicated, one may show that the 
conclusion remains unchanged. Thus there are no solutions of the BPS equations for Ov 2 > 1 
that can be expanded in a power series of u> n . Here we do not turn on the diagonal entry A a 
of K. As will be explained later, the effect of nonzero A a corresponds to locating each vortex 
at A a = A" — i\ a y position. Considering the case of one vortex, one can easily turn off this 
value by using the translation symmetry of the system. Hence our proof above is strictly 
applicable to this case. Furthermore, m vortices are an assembly of individual vortices, one 
naturally expects that the above proof goes through even m vortices with generic values of 
A,, 

Figure 1 summarizes our investigation of the static solutions in the Abelian Higgs 
theory for A = 1. The self-dual BPS solutions exist only for Ov 2 < 1. In the range the non- 
BPS exact solutions are unstable due to their higher energies. When Ov 2 > 1, the non-BPS 



branch alone continues to exist. For the solutions of a negative flux, it is shown in Ref. [10 



that the anti-self-dual solutions exist for Ov 2 ^> 1 or Ov 2 <C 1. In the intermediate values of 
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9v 2 , the existence of the self-dual solutions is not known. The exact vortex solutions with 
positive magnetic flux exist even for A ^ 1. It will be shown later that they are also stable 
only when 9v 2 > 1. 



4 Fluctuation spectra around the vortices 



In the last section, we have identified the possible static vortex solutions including BPS 
and non-BPS cases. In the BPS case, the classical stability of the solution is quite clear 
because the energy is saturating the bound set by the topological quantity. For the case 
of non-BPS, however, it is not a priori clear whether the vortices are stable or not. When 
9v 2 < 1, we have shown that there exist solutions that have lower energies than the exact 
non-BPS solutions. Thus we expect naturally that there should be tachyonic modes. It 
is also shown that BPS solutions do not exist for 9v 2 > 1. Hence in this case the issue 
of stability seems a different matter. To resolve these issues clearly, we shall study, in this 
section, the quadratic fluctuation spectra around the exact solutions identifying the signature 
of mass squared for all possible degrees. It turns out that the solutions 9v 2 < 1 is indeed 
unstable by developing tachyonic modes in their spectra. In case of 9v 2 = 1, the potential 
tachyonic degrees become massless and the solution is indeed stable. For 9v 2 > 1, solutions 
are classically stable because the tachyonic degrees become massive. For all these three 
cases, the first m diagonal elements of the gauge field fluctuation are massless, which will be 
identified with the degrees of vortex positions. 

Let us study first the quadratic fluctuation of the original theory about the vacuum 
K = c + /C and = v (1 + h) without any vortices. The Lagrangian is then reduced to 



2tt 



-tr 



-2v' 



\]C\ 2 + 9v 2 (\h R \ 2 + \h I \ 2 )-±\[c,tf] + [lC,J 



29 



-2v' 



[c, h R ] 



K, + i[c, hi] 



2\9v 4 h 2 R 



with the gauss law constraint, 

[c\K] + [c, /C f ] - 2i8v 2 h I = , 



(34) 



(35) 



where h R = + frf) and hi = ^(h — ftf). One may simplify this action by reintroducing 
A Q field, which has a role of imposing the Gauss law constraint. We then choose a gauge 
A — hi, at which /C + i[c, hi] — > K. The Lagrangian becomes 



T 

L = — rtr 



JC\ + 9v 2 \h R \ 2 - Y e \[cX] + [/C,ct]| 2 - 2v 2 (\[c, h R ]\ 2 + |/C| 2 + X9v 2 h 2 R ) 



with the gauge condition now 

[c\ t] + [c, t ] ] - i[c\ [c, hi]] - i[c, \c\ hi]] - 2i9v 2 hi = . 



, (36) 



(37) 
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This can be solved in terms of hj for arbitrary /C, on which the Lagrangian does not depend. 
It is now clear that all the degrees are massive; the components of /C have a mass squared 
greater than 2v 2 while Hr components a mass squared greater than \v 2 . We see that the 
gauge field absorbs part of the scalar degrees and becomes massive. This corresponds to the 
so called Higgs mechanism of the ordinary gauge theory when the gauge symmetry is broken 
spontaneously. 

To study the quadratic fluctuation around the exact solutions, we turn on generic 
perturbation of the form 



K = SmcSb + A + K , <f> = v(l + (p)S m , 
with /C and if decomposed as 



(38) 



JC = A + VSl + s m w^ + s m icsl 

( o X \ 



A V 

it 



(39) 



<P = XS* m + S m hSl 



h 



Here we set (pP m = with out loss of generality since an arbitrary 5<j) can be expressed by 
vipS m . Further introducing unitary operators 



^(A a c- A„ct) 



U a = e 



we parameterize the components deliberately as 

m— 1 m—1 



(0 < a < m - I) , 



(40) 



m—l oo 



A =J2J2 A ab\a)(b\, 

a=0 6=0 

m—l oo 

W = Y,H W an\a){n\U a , 



V = J2J2 V an\a)(n\U a , 

a=0 n=0 
m—l oo 

X = Y,H X an\a){n\U a , 



(41) 



a=0 n=0 

oo oo 



a=0 n=0 

oo oo 



£ = J2^kn \k)(n\, 

k=0 n=0 

The unitary operators, U a , satisfy 



h =J2 J2 h kn I k ) ( n | . 

k=0 n=0 



UJU\ = c t + -^A a , 



(42) 



which are helpful in identifying variables that diagonalize both the kinetic and potential 
terms. 



Now we insert these to the original Lagrangian in the gauge A = and expand it to 
the quadratic terms of the fluctuation. We get 



T -- 

9 2 

q 2 ^ 

if a,n 



^(|Co6| + \Qab\ 2 - -^-^-^ \Cab\ 2 ) + ^2[\T a \ 2 - 



ab 



\H n 



+ \Y an \ 2 + \G 



e 2 

2 2n + l + 9v 2 
6 



9v 2 -l 
6 



(\H I 2 + \Y \ 2 ) 

l I- 1 - 1 an | ~ | 1 an | J 



T a 

+ L D , (43) 
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where Ld is same as the Lagrangian ([34]) but /C and h and are replaced respectively by /C 
and h. In this Lagrangian, we put 

C ab = ± [e-^A ab - e ie ^A[ h ] , Q ah = ^ [e~^ At + e*-^] , (44) 

where a (, is the argument of A a — A;, or an arbitrary constant for A a = A&. We also set T a , 
Han-, Y an and G an as 



Vth 2 



H a n — rj. — y, 

y/6v 2 + 2n + 



= [X an y/2n+l - J— (Vn+1 V a , n+l + VnW a ,n-l)) , 

1 \ y/2n + l v y / 



V2n + 1 v 



1 



G an = .„ = = = (^V 2 X an + V^+T K, »fl + V 7 ^ tVo, r^l ■ (45) 

To this order, the Gauss law constraints for A and the off diagonal degrees become 

g ab = (only for A a ^ X b ), G an = , (46) 

and, for Ld, it takes the same form in fl3"5]) where /C and /i are again replaced respectively 
by JC and /i. 

From this it is clear that C a b is massless when A a and X b coincide. In particular the 
diagonal components C aa and Q aa are always massless; they are associated with the transla- 
tional motion of the vortices. The nature of this motion will be exploited when we discuss 
the low energy dynamics of the vortices. When 6v 2 < 1, T a has a negative mass squared. 
Hence we see that the vortices are unstable even for the case of a vortex. On the other hand, 
for 9v 2 > 1, the instability disappears and the vortex solutions are stable. This is also quite 
consistent with the fact that there are no BPS solutions for 9v 2 > 1. If there were such 
solutions, there must be tachyonic modes because the BPS solution should have the lower 
energy than the non-BPS solutions. 

Especially when 9v 2 = 1, the potential tachyonic degrees become massless and may 
participate in the low energy dynamics as will be discussed later. The remaining of diagonal 
components are H an and Y an . The G an degrees are dropped out of the physical space spec- 
trum once the Gauss law constraint is imposed. Here we were be able to diagonalize these 
infinite dimensional degrees, which is in general not an easy task to achieve. The spectrum 
of these physical degrees is particularly simple; they are all massive with the same mass 
v 2 + ^j^i which is independent of the index a. This spectrum can surely be understood 
from the underlying D-brane perspective. 

Finally, Ld describes the fluctuation spectra of the original system around its trivial 
vacuum configuration. This is no coincidence because the degrees of the original system still 
remain around vortices. At this point, we like to emphasize again that they are all massive 
controlled by the mass scale v and y/Xv . 



10 



5 Low energy dynamics 



From the analysis of the fluctuation spectra, it is clear that the vortices are unstable due to 
the tachyonic modes for 9v 2 < 1. On the other hand, the vortices do not exhibit any tachyonic 
instabilities for 9v 2 > 1. For all ranges of parameter 9v 2 , the vortex solutions depend upon 
2m-dimensional free parameters where m is the topological number corresponds to the total 
number of vortices. We shall first consider the stable case where 9v 2 > 1 and begin by 
clarifying the physical interpretation of these parameters. In short, these parameters, A a , 
are positions of vortices on two plane where the noncommutative gauge theory is defined. 
For the gauge group element defined by 

U P U P = I (47) 

with P m Up = UpP m = P m , the corresponding gauge transformation affects only first m x m 
and m x oo component of K and <fi. Utilizing this gauge freedom, we have diagonalized 
m x m part of K by 

P m KP m = —7= diag[A , Ai, • • • , A m _i] . (48) 
V 29 

in the solution (0). Any permutations of the eigenvalues A a and A& are achieved through 
the gauge transformation by the Weyl subgroup elements. So they are physically equivalent 
configurations. Thus the moduli space is in fact (R 2 ) m /S m where S m is the permutation 
group. 

In order to identify the meaning of the moduli parameters, let us first study the ef- 
fect caused by the overall translation of the vortex solutions. For this, we note that the 
infinitesimal translation is given by 

= - i^A^) = -bDrf , (49) 

where we have added the infinitesimal gauge transformation by the gauge function £jAj. On 
the solution, this produces 

SAi = ^ei&Pm, 50 = . (50) 

The magnetic field B and the Higgs gradient are unchanged by the translation and, 
consequently, one may construct easily the fields translated by a finite amount. Namely, 
the Higgs change is A0 = while the change of gauge field in terms of K variable is given 
by AK = -j=£P m where £ = £i — i^- Hence we see here that the total translation leads 
to a uniform shift of each A a by the amount £. This is of course quite consistent with the 
interpretation that the moduli parameters represent positions of the vortices. Of course due 
to the U(oo) gauge symmetry, the effect of translation does not quite look like a translation 
of profile in case of ordinary field theory where a density p(x), for example, is merely shifted 
by £ as in p(x — £) as a result of the translation. In this respect whether local informations 
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such as positions of vortices is well defined in noncommutative gauge theory is not obvious at 
first sight. There is another way to get the above result of translation. The global translation 
generator can be alternatively expressed as 

T = e~^ n . (51) 

where Pi is the translation generator p$ = — ^[e^Xj, • ]. In noncommutative field theory, the 
operation of translation on a field can be expressed as a similarity transformation, 

Tf(x) = U T f(x)tf T . (52) 

where Ut is a unitary matrix defined by ev^^ c t ~^ c \ In our case, we add a gauge transfor- 
mation by U = Ut after the translation. Then the resulting gauge and scalar fields read 

A' = A + i{, 

0' = U T • (53) 



The gauge field is shifted only by a constant piece. We see also that (jxjr is invariant. If the 
scalar were in the adjoint representation, it would be invariant under the transformation. In 
order to obtain the previous result in (pUD, we further perform a gauge transformation by 

U = e^^ 1 ™ ^ Cm ' ) with c m = S m cS m . Using the explicit expression of the solutions, one may 
easily check that results agree with 



One could also study the exact solutions moving in a constant velocity as discussed in 
Ref. [0. The theory is not Lorentz invariant because the ^-product does not respect the 
Lorentz symmetry. However, as discussed in Ref. JF|, one may still construct moving soliton 
solutions once the static solution is given. The construction is achieved by Lorentz boosting 
of the static solution followed by the change of 9 by 7 where 7 is the Lorentz dilation 
factor defined by l/\/l — f3 2 with velocity (3. Constructed this way, the solution moving in 
x direction reads explicitly, 

A o = -lf3xA x (x', y'\ 7 0) , A' x = ^A x (x', y'\ 76) 

A' y = A y (x', y'; 7 0) , <t>' = <f>(x', y'; 7 0) , (54) 

assuming Aq = for the static solution. Here the arguments are given by 1' = 7(1 - j3 x t) 
and y' = y and the fields without prime denote any static solutions. In the present case, 
one may further simplify the form of the moving solution again taking A' Q = 0. This gauge 
choice is achieved from the above solution by the gauge transformation with 



where we define 



u = e -^^si-is m csi) e ^(z#-m^ (55) 



'20 

00 

S = J2 \n + = c^cct)- 1 , (56) 



n=0 
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and £ = j3 x t. Here \n)' is the number eigenstate constructed by the number operator c'c. 
The form of the solution becomes^] 

^4 = —i—j=(c — SmCSjn) + — j3 x t P m 

<P' = vS m , (57) 

with A' Q = 0. The map from c = ^= to the new basis c = h= (x — iy^yfl)*) belongs to 
the area preserving diffeomorphism. Except some overall numerical coefficients, the solution 
apparently represents a configuration that has an elliptic shape; for example, the magnetic 
field of a moving vortex appears in the function representation as 2e~e( x f+v t ). Utilizing 
the U(oo) gauge symmetry, the solution Q57p can be further mapped to 

^4 = —i—j={c — SmcSjn) + —/3 x t P m 

<P' = vS m , (58) 
by the gauge transformation with the unitary matrix 

m—1 

U s = S m Sl + £ \a)'(a\. (59) 

a=0 

Inserting (|58|) into the time dependent field equations, one may directly check that it is 
indeed a solution. Actually, one may even construct solutions representing more general 
motion of vortices. The time dependent solutions read 



m—1 



A' = -i^=(c - S n cSl) + - £ i xa + P^] l«><«l 

</>' = vS m , (60) 

with (3 a = /3™ — The motion of each vortex takes place independently to an arbitrary 
direction. The magnetic field and electric field are now 

i m—1 

B' = -P m , E , i = e ij y £ft\a)(a\' (61) 



3 

a=0 



The energy of the moving vortices is evaluated as 

•m-SHsw-^-^), (««) 

where no approximation is made. The energy behaves precisely as free nonrelativistic parti- 
cles with a mass 



^The appearance of c instead of c in the gauge field is not a typographical mistake. 
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One striking fact in the moving solutions lies in the fact that there seems no limit in 
the velocity. It can apparently exceed the light velocity!]. On the other hand, in the original 
construction by the Lorentz boost followed by the change in the scale 9, the construction 
itself loses its validity when the velocity exceeds the light velocity. Specifically, the factor 7 
becomes imaginary. Nonetheless, the final form of the solution in this range of velocity does 
solve the time dependent equations of motion. Our system lacks the Lorenz invariance and, 
thus, this seems not a serious trouble. But one should remember that the system originally 
has the Lorentz invariance while a specific background field (the constant NS-NS two form 
background field) is turned on in the string theory context |L9| . Without going into details, 
we like to mention the fact that, when the velocity exceeds the light velocity, part of once 
stable degrees become tachyonic and instabilities are necessarily set in. Hence the solutions 
seems not to have much physical significance when the velocity exceeds the light velocity. 
Further investigation is required on this issue. 

Let us now turn to the moduli dynamics of vortices. The study of translation justi- 
fies that A's faithfully represent overall position of vortices. Let us consider the following 
operator, 

Xi = Xi - OtijAj , (63) 

which may be rewritten equivalently as X = Xi — iX 2 = vWK. This transforms covariantly 
under the gauge transformation, i.e. X — > WXU. Since the operator reduces to Xi in the 
commutative limit and is gauge covariant, we shall call it as covariant position operator. 
Another justification for the terminology comes as follows. It transforms as 

X[ = X i + i (64) 

under the translation of (|52[) followed by the gauge transformation by U = Ut- This is 
precisely the required property as a position operator under translation up to gauge freedom. 
It will be used to measure local properties of the noncommutative field theory. To show that 
the eigenvalues A a represent positions of vortices, let us consider the following moments 

I ki i = 27r6tr[X k (X r ) l H] . (65) 

These quantities are gauge invariant and measure the local distribution of matters in non- 
commutative gauge theory. For example, ii.i corresponds to the moment of inertia for the 
configurations of the ordinary field theory. 

For the exact vortex solutions (6v 2 > 1), we have 

m— 1 

hi = M one £ \ k a \ l a . (66) 

a=0 

II If the moving solution were not exact, we would have easily missed this point. This is similar to the case 
of the noncommutative scalar field theory with a quartic interaction. The two particle bound state energy 
is unbounded from below, which was observed in the exact nonperturbative computation of the bound state 
energy |J|. 
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In case of commutative field theory limit, the same moments can be found only when the 
Hamiltonian density is sum of delta function as H(x) = J2T=o M one 5 2 (]t ~ A ). Thus we 
show that the relatively local information of noncommutative gauge theory can be obtained 
from the moments defined above and that the eigenvalues A a are representing the positions 
of vortices up to the permutation symmetry. Considering, for example, vortices located at 
the origin, the size information of the vortex configuration can be extracted by the moment 
of inertia. The "size" (measured by the covariant position operator) is finite for the BPS 
vortices (6v 2 < 1). In fact it decreases within the BPS branch as 9 gets larger and becomes 
zero for the stable non BPS vortices (9v 2 > 1). 

The moduli dynamics of the noncommutative solitons may be pursued in a similar 
manner as solitons in a ordinary field theory. As stated before, we shall consider first the 
case where Ov 2 > 1. We proceed by giving the time dependence to the moduli parameters and 
adding an appropriate gauge freedom so that the motion respects the Gauss law constraint. 
But in our present case, it is enough to simply give the time dependence without adding any 
gauge degrees because they already satisfy the Gauss law constraint. Namely we insert 

K = K(X a (t)) , = 4>(\ a (t)) (67) 

to the full Lagrangian where quantities with bar denote the vortex solutions. (This ansatz 
is quite consistent with the moving solutions constructed before.) The resulting effective 
Lagrangian is given by 

71 m_1 • - 

L eS = -mM onc + — • (68) 

Consequently, the moduli space metric on (R 2 ) m / S m is flat, i.e. 

TO— 1 

ds 2 =Y / d ^d\ a . (69) 

a=0 

The inertia mass here is different from the rest mass but there is no physical reason why 
these two masses agree. Not to mention, this effective Lagrangian can be easily quantized 
and wave functions are those of m free nonrelativistic bosons with mass 4n;- 

9 8 

For the present model, one may in fact go beyond the moduli space description in dis- 
cussing the relevant low energy dynamics. In the previous section, we studied full fluctuation 
spectra around the static solutions. We find that off diagonal degrees W a j, V a j (j ^ 0) and 
P m (pP m (with 4> = v (1 + (p)S m ) are massive with a mass squared m 2 = ^±1 + v 2 (j > 0). 
Furthermore P m lCP m components have a mass squared at least order of v 2 . The real part 
of P m if P rn has a mass order of Xv 2 while its imaginary part is a gauge degree of freedom 
that will be absorbed into the gauge field P m K,P m by the Higgs mechanism. The alternative 
description of low energy dynamics is obtained by ignoring all these massive degrees of free- 
dom and focusing on all the remaining fluctuations around the A a = solution. Namely we 
only consider the fluctuation of the gauge field in m x m sector and the potential tachyonic 
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mode defined by 



A = P KP 

m—1 



(70) 



a=0 



The full Lagrangian is then reduced to 

To 



Left = 



eti AA ] - hi\A,A ] } 2 + e\ \t) 2 1 



•A t |r>| 2 + |.Ak>| 2 ) 



-^(r\[A,A\r)-(ev 2 -l) 





2 




4" 


\r) 




\r) 





(71) 



where the Gauss law constraint 



[A, A 



- [AA ] } + 
A ] \f) = 



\t)(t\ - \t)(t\ = 



(72) 



is still in effect on the Lagrangian. The residual part of the Lagrangian can be organized 
as follows. Denoting all the remaining massive modes collectively by Z p , there are terms of 
O(Zp), O(AZ^), 0(tZ%), 0(tAZ p ), 0(t 2 Z p ), O(Z^) and quartic terms including at least 
one massive degrees Z p . One should note that there are no terms of order 0(A 2 Z p ). 



When 9v 2 > 1, the tachyonic modes become massive too. To truncate the Lagrangian 
consistently, we consider A ~ 0(e). The the A 4 terms contribute to the Lagrangian as 
0(e 4 ). Now if one turns on any massive degrees, it should be 0(e 2 ) due to Z p or r 2 terms in 
order to have a valid approximation of dropping the massive degrees. Then the interaction 
terms between the massive and the massless degrees are of higher order, i.e. 0(e n ) with 
n > 5. For example, we see that the terms of 0(AZ 2 ) is of order 0(e 5 ). If there were terms 
of order 0(A 2 Z P ), these would contribute to the potential as 0(e 4 ). However, there are no 
such terms as stated previously. Hence the massive degrees are effectively decoupled from 
the massless degrees to the quartic order in the low energies. Hence we may consistently 
drop all the massive degrees consistently. Ignoring all the massive modes, we are led to 



J eff 



27T 

7 



2o 



with a constraint, 



[A,A t ] + [A\A]=0. 



(73) 



(74) 



This Lagrangian is precisely the matrix model, which coincides with the bosonic part of an 
effective Lagrangian for m DO-branes moving in two dimensional target space. 

The vacuum moduli of this effective action is the vortex moduli described previously 
by the coordinate A a on (R 2 ) m jS m . We see clearly that the singularity when vortices are 
overlapping is resolved in this description. Moreover, the commutative moduli coordinates 
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are replaced noncommutative matrix degrees whose structure is especially relevant when 
vortices are nearly coincident. Hence legitimate approach toward the quantization of the 
low energy dynamics is also quite clear. 

One might ask at this point about the nature of the coordinates of vortex positions. 
Since the noncommutative space underlies in defining the noncommutative field theories, one 
would also expect that the noncommutative solitons should see directly the noncommutative 
nature of the underling space through their forms of interactions. But above description does 
not show directly the noncommutative nature. Namely, the interactions does not show any 
particular structure depending upon the noncommutativity scale 9. Stated again, nothing 
particular happens at the separation AA ~ \J~Q~. Nonetheless, the vortex positions are truly 
described not by c- number eigenvalues but by matrices. In this respect, the locations of 
vortices still possess a noncommutative nature that is originated from the matrix properties. 

Next, we consider the case where 9v 2 = 1. In this case the potential tachyonic modes 
become massless. But there exist a quartic contributions, so it is not a moduli degree 
as defined by the configuration space of the constant energy. But we include it because 
its contributions is of the same order of A when A is small. Hence, to study interaction 
between the massless modes and the massive modes, we let A and r be order of e as before. 
Then Z p may be allowed to the order of e 2 to have a well defined low energy description. 
But this time, there are interaction terms of the forms tAZ p and T 2 Z p , whose contribution 
to the Lagrangian is 0(e 4 ). Hence the massless degrees are not decoupling from massive 
degrees. One could write down the consistent effective Lagrangian for this case too. But it 
turns out that the effective Lagrangian involves infinite number of massive degrees. Instead 
of giving detailed analysis, we here briefly comment on the nature of the resulting motion 
involving the potential tachyonic modes. First note that one may effectively describe the 
motion by A and r once all the massive modes are integrated out. One may then easily 
verify that, among 0(A 2 t 2 ) terms, only the term of (r | [A, A'] |r) remains out of (|7T|). This 
is quite consistent with the translational invariance of the underlying system, whose action 
is replacing A by A + £/ mX m- The motion along r-direction is controlled by two terms, 
(t|[/4., *4J]|t) and r 4 . This motion excites other components of the magnetic field out of the 
static solution |P m , while the flux $ preserved. Hence the motion represents an oscillatory 
dispersion of magnetic field to other components. If the tachyonic modes are small enough, 
the part of matrix mechanics responsible for the vortex positions is, is little affected for fixed 
energies. 

Now we turn to the case where 9v 2 < 1. In this case the fluctuations include the 
tachyonic modes. Small fluctuation will trigger the vortex to run into a more stable lower 
energy configuration that corresponds to BPS states. As shown previously, the BPS state 
has the same flux as the original unstable static configuration. Thus during the process, the 
flux should be conserved while the difference in energy is eventually dissipated away. The 
tachyonic instability is present even for the case of single vortex. So it can be interpreted as 
a collapse of each individual vortex to a more stable one, i.e. the BPS state. The detailed 
study of the collapse will be quite interesting in relation with recent discussion of the tachyon 
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condensation in string theory. 



6 Conclusions 



In this note, we have first investigated general static soliton solutions in the noncommutative 
Abelian-Higgs theory. There are exact multi-vortex solutions found in Ref. |J for general 
values of parameters A and 9v 2 . These are in general non-BPS except A = 9v 2 = 1. We 
extend these solutions by finding exact solutions describing vortices positioned at arbitrary 
locations. We have shown that these solutions are unstable only when 9v 2 < 1. It is therefore 
expected that lower energy non-BPS solutions exist for 9v 2 = 1 and A ^ 1. We confirm this 
by considering self-dual BPS branch for A = 1. For 9v 2 < 1, the self-dual BPS branch 
develops, which has an lower energy than the exact unstable vortices. The BPS branch 
ended at the point 9v 2 and there no longer exist BPS solutions for 9v 2 > 1. Instead, the 
exact non-BPS configurations become stable configurations. We also illustrated the case of 



anti-self-dual BPS solutions that have an negative flux|10||. The solutions are shown to exist 
for 9v 2 <C 1 or 9v 2 ^> 1. For the intermediate region, the existence of the BPS solutions are 
not clear yet. 

We then discussed the general fluctuation spectra around the exact static vortices with 
general moduli parameters A a . It is shown that there are tachyonic instabilities only when 
9v 2 < 1. We have identified the massless degrees of freedom and masses of all the off diagonal 
degrees. With help of the covariant position operator and studying translation of vortices, 
we were able to identify the physical meaning of the moduli parameters; they are positions 
of the vortices. We were able to construct exact moving solutions of vortices, where each 
vortex is moving freely in a arbitrary constant velocity. We then show that the metric in 
the moduli space is indeed flat by evaluation the low energy effective Lagrangian within the 
moduli space description. In fact one may go beyond the moduli space description in this 
case by identifying quartic order interaction terms of the massless degrees of freedom. It is 
nothing but the matrix model of m DO-branes moving in a two dimensional target space. 
Thus we have shown that the low energy dynamics are faithfully described not by positions 
of individual vortices but by matrices. 

The exact time dependent solutions describe vortices with constant velocity. What is 
striking in the solution is not that vortices are moving freely but the velocity is not limited by 
the light velocity. The solution exists even for the velocity greater than the light velocity. We 
argued that the fluctuation becomes tachyonic when the velocity exceeds the light velocity. 
Therefore, the solution seems not to have much physical meaning when the velocity exceeds 
the light velocity. Further detailed study is required on whether the solution in the region 
is consistent with special relativity or not. Though the system lacks the Lorentz symmetry, 
the special relativity should be still in effect because one may regard the system as Lorentz 
invariant system with a specific background field (a constant NS-NS two form background 
field in string theory) is turned on. 
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We expect that our investigations can be generalized to the N=2 supersymmetric ver- 
sion of the noncommutative Abelian-Higgs theory. In particular, supersymmetries will not 
be preserved even partially for the sector of nonvanishing flux with 9v 2 > 1. We like to 
finally mention that our investigations may be applicable to other exact solutions recently 
found@,|T3|,|T§,|r7|. 
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A Self-dual BPS solutions 



Here we demonstrate how to obtain the BPS solutions. Substituting the ansatz for K 
into the master equation (|^) gives 



J2(ev 2 -l + \f n \ 2 -\f r ^\ 2 ) \n)(n\ 



n=0 

oo oo oo 



EEEt^ f**"' fnin-fyfj f#P'" fiM.n-i) P \i+np) (j+np\ , (A.l) 

i=0 j=0 n=0 n - 



where we set \4>q) = X^o s «K) an d fj = f° r an y 3 < 0- Comparing | 0)(i |, i > 1 
components of the left and right sides we see = soSj, i > 1. Now by mathematical 
induction, one can show easily = SiSj, i ^ j. Hence we may put \ (f> ) = ( \ m) ioi some 



complex number ( and a non-negative integer, m. This simplifies Eq. ( A.l 



as 



£(V-l + |/ n | 2 -|/ n g 2 ) \n)(n 



n=0 

oo 



\C\ 2 J2 —\ I/™! 2 IfrrMvl 2 - ■ ■ \fmM,n-i)p\ 2 \m+np)(m+np\ . (A.2) 



n=0 n! 



Hence 



9v 2 -l+\f n \ 2 -\f n ^ p \ 2 =0, (A.3) 
for < n< m or m < n, n ^ m mod p, and 

@ v2 — 1 + \ frrHnp\ 2 ~ | /m+{n-l)p | 2 = ^rl^ m | 2 l/™+p| 2 ' ' ' \ frn+{n~l)p\ 2 ■ (A-4) 

for < n. With q n = l/^^pl 2 — n — 6v 2 {p — l)(n — 1), the magnetic field is expressed as 

oo 

B = v 2 ^2' \n)(n\ + [Qn - q n +i - Qv 2 {p - 1)] | np + m)(np + m \ , (A. 5) 

n " n=0 
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where J2' is the sum over < n < m and m < n, n ^ mmodp. The flux is then given by 

$ = m 9v 2 + go — lim q n , (A. 6) 

n^oo 

where if we write m = pk + r, < k,0 < r < p — I, q = k(l — 6v 2 ) + (p — l)6v 2 . In order to 
have a finite energy q n ought to converge. If ( = 0, we find that \f np+r \ 2 = (1 — 0v 2 ){n + 1), 
< n, < r < p and q n = —9v 2 pn + 6v 2 (p — k — 1) + k. For this solution, the energy 
diverges. On the other hand, if £ ^ then q n satisfies the following recurrence relation 

~ 7T~2 = 1 + (P - + - 9n - (P - 1 ^ • A.7 

We take the n — ► oo limit of the above equation and conclude that p = 1 is a necessary 
condition for g n to converge. Now for p = 1 let us assume that limg n = a. This implies 
that, for any e > 0, there exists large N such that a — e < q n < a + e for n > N. Eq. ( |A.7| ) 
implies 

— 9jL -f^ = n (i + -) • < A -s) 



Furthermore, we have 



n(i + ^)<n(i + |)<n(i+^). ^ 

n=iV " n=iV U n=N n 

However for any g ^ 

oo / oo \ 

H (l + -J = exp H n + e) -Inn , (A.10) 



n=iV ' n \n=7V 



which is either infinity or zero depending on the signature of g. Thus, Eqs. ( A.8 - A..SQ implies 
that a must be zero. 



With p — 1, Eq. (|A.7 ) gives a recurrece relation 

q n +i -q n = q n - Qn-i + —(<?«- Qn-i + 0v 2 ) , (All} 



n 



with two initial data, go — m(l — 6v 2 ) and q\ = m(l — 6v 2 ) + \(\ 2 — 8v 2 . Therefore if \(\ 2 > 6v 2 , 
then q n is monotonically increasing. As the only possible converging value is zero, it must 
diverge. In case \(\ = 0, it can be easily solved by q n = m(l — 6v 2 ) — n6v 2 . 



B Nonexistence of self-dual BPS solutions for 6v 2 > 1 

We shall work in a gauge = for i > j. The BPS equations can be written as 

(pel = h ] P m + c ] m ipP m + h ] V P m (A.l) 
tP m = (l + e)(<pP m + P m ^ + <pP m ^)-([c m ,tf] + [h,cl} + [h,tf})^ (A.2) 
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where c m = SmcS^ and P m = 1 — P m . The relevant part of the first order equations in u 
reads 

Pm { fi(\)Pm = Pmh(i)Cm PrrJ^{X) < ^rn 

PmV{l)4n = Pmh\x)P m • (A. 3) 

Since P m h^P m = for our gauge choice, we find that = h^ m+1 , foW can be arbitrary 
but all the remaining components should vanish. Now we investigate the second order 
equations obtained from the perturbation equation ( |A.2p . Let us multiply P m to the left 
and to the right of the equation at the same time. We obtain 

— n-Pm + Pm^(l)Pm i p i [i^Pm = \Pmh(l)P m , Pm.h(l)Pm] + Pmh(l)P m h^P m ■ (A-4) 

Using the result of the first order equations and taking trace of the above equation, one finds 

me m ~ x 



-TT=EI^| 2 - (A.5) 



a=0 



Hence we get a contradiction when e > 0. 
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